# Ref: from NIST website

from Excitation_Hydrogen import *

cs = CS_Excitation(species_name = "H")

for id in np.arange(2,11):
    reaction_name = "Excitation_H-1s-" + str(id) + "p"
    cs.read_cross_section(reaction_name)

cs.plot_cs()
cs.save_cs()
cs.cal_total_cs()
